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Explicit analytical expressions for the drag and diffusion coefficients of a spherical particle 
attached to the interface between two immiscible fluids are constructed for the case of a 
small viscosity ratio between the fluid phases. The model is designed to explicitly account 
for the dependence on the contact angle between the two fluids and the solid surface. The 
Lorentz reciprocal theorem is applied in the context of a geometric perturbation approach, 
which is based on the deviation of the contact angle from a 90°-value. By testing the 
model against experimental and numerical data from the literature, good agreement is 
found within the entire range of contact angles below 90°. As an advantage of the method 
reported, the drag and diffusion coefficients can be calculated up to second order in the 
perturbation parameter, while it is sufficient to know the velocity and pressure fields only 
up to first order. Extensions to other particle shapes with known velocity and pressure 
fields are straightforward. 


1. Introduction 

The diffusive behavior of colloidal particles is drastically altered compared to diffusion 
in a bulk fluid when the particles are affected by the presence of an interface between two 
immiscible fluids. The motion of particles attached to a fluid interface occurs predominantly 
parallel to the interface, but may also involve temporary particle detachment, as can be 
concluded from experiments ( Walder et al.|[2M0 Sriram et q 1.||201'2 |. The phenomenon 
of two-dimensional interfacial diffusion is not yet fully understood, which is reflected in 
the variety of experimental results as well as related theoretical models. For example. 


Peng et al. (2009) and Chen & Tong (2008) have studied the influence of the surface 


concentration of diffusing particles on the diffusion coefficient. In the limit of infinite 
dilution, the measured diffusion coefficient is found to be very close to the bulk value in 
one of the fluid phases. The authors of both publications explain the data by assuming 
the interface as incompressible, although there is no reason to assume contamination of 
the interface ( Peng et aL||2009 ). As another example, experimentalists have occasionally 
found the size depen dence of the diffusion coefficient t o differ from the inverse of the 
particle radius, a~^ (Du et al. 2012 Wang et al. 2011), at variance with the modified 
Stokes-Einstein relation (Brenner fc Leal||1978) 


D = 


kT 


67r^ia/(0,I^2/Mi) 


( 1 . 1 ) 


By the index 1 we denote the fluid with the higher viscosity /ii, while index 2 denotes the 
fluid having the lower viscosity. The three-phase contact angle 0 is measured in fluid 1. 


In equation (1.1), the function / specifies the deviation of the drag force from the Stokes 


drag of a spherical particle suspended in the bulk of fluid 1. In terms of the drag force Fjj 
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acting on the attached particle, / is thus defined by 

Fd = -Gnti^afU, ( 1 . 2 ) 

where U denotes the particle velocity relative to the undisturbed fluids. Even when the 


modified Stokes-Einstein relation (1.1) is valid, as we shall assume in this study, the 


functional form of the drag coefficient / of a translating interfacial particle is not known. 
A variety of theoretical and experimental studies deal with the drag coefficient of particles 


attached to fluid-fluid interfaces (Fulford & Blake||I986 O’Neill et al. 1986 Petkov et al. 

1995||Danov et a/.||I995 

1998 Cichocki et aL||2004| Fischer et a/.||2006 Pozrikidis||2007| 

Ally & Amirfazli|2010 I 

Ilawzdziewicz et aZ.||2010 

). However, the corresponding theoretical 


models mostly rely on numerical methods. With this work, we intend to contribute to the 
field by providing an explicit analytical expression for the drag coefficient of a spherical 
particle attached to a pure interface between two fluids of highly different viscosity. 
According to the Stokes-Einstein relation (O), the diffusion coefficient D directly follows 
from the drag coefficient /. Our analysis aims at the latter, which can be studied by 
means of low Reynolds number hydrodynamics. The main feature of the new model is 
to account for the dependence of the drag and diffusion coefficients on the three-phase 
contact angle. 


2. Series Expansion of the Flow Field 

Our modeling focuses on the drag coefficient of a rigid sphere translating along a fluid- 
fluid interface at low Reynolds number. We hereby study the fundamental case of a pure 
fluid-fluid interface. Therefore, we do not employ any incompressibility constraint for the 
interface and assume a vanishing interfacial viscosity. Also, we neglect any deformation of 
the fluid-fluid interface on the scale of the particle radius a, corresponding to a negligible 
influence of external forces acting normal to the fluid-fluid interface, such as buoyancy 
or electromagnetic forces. Assuming a planar interface also implies that the capillary 
number Ca := fiiU/a (with U = ||J!7|| and the fluid-fluid interfacial tension a), measuring 


dynamic deformations of the fluid-fluid interface, is small compared to unity (Radoev 


et al. 1992). In addition to interfacial deformations, we neglect particle rotation. The 


validity of this assumption depends on the conditions at the three-phase contact line, 
for which two limiting cases exist. Firstly, the contact line can be pinned to defects of 
the particle surface, meaning that it retains a fixed position with respect to the latter. 
A pinned contact line prevents the particle from rotating if the capillary forces due 
to the fluid-fluid interface are sufficiently large. Secondly, the contact line can move 
tangentially to the particle surface if the contact-angle hysteresis is small or vanishes. In 
this case, the particle may rotate with an angular velocity dependent on the rate of 
dissipation occurring at the contact line. If the angular velocity is negligible compared 
to U/a, the particle may be approximately considered non-rotating, because then the 
particle’s surface velocity associated with the rotational motion is much smaller than 
the surface velocity U associated with the translational motion. For increasing angular 
velocity, that is, decreasing dissipation at the moving contact line, the assumption of a 
non-rotating particle looses its validity. Consequently, the following model is valid for a 
particle with Da/U ^ 1. To further simplify the mathematical treatment, we assume a 
vanishing viscosity ratio between the two fluids, /i 2 /Aii 0. Thus the planar fluid-fluid 
interface effectively becomes a symmetry plane from the viewpoint of fluid 1, while the 
influence of fluid 2 can be neglected. For this reason, the flow problem is analogous to the 
motion of a body possessing reflection symmetry moving parallel to its symmetry plane. 
Clearly, for a contact angle 0 of 90°, the symmetric body is spherical, resulting in the 
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classical Stokes flow problem around a sphere. The corresponding drag on an interfacial 


particle is then simply half the Stokes drag in the bulk of fluid 1 (Radoev et al. 1992 
|Danov et a/.||1995 Petkov et al.||199^ Ally fc Amirfazli]|2010 ), implying 

/(90°,0) = l/2, (2.1) 

where the function / is defined by equation For contact angles differing from 90°, the 
symmetric body consists of two fused spheres (cf. figure [^b)). This case has been studied 
by Zabarankin (2007), who provides drag coefficient values derived from a numerical 
solution of a Fredholm integral equation. Recently, |D5rr fc Hardt (2014) have derived the 
asymptotic expression 


/(0,O) = - 


1 H-cos 0 + C>(cos^ 0) 

16 


( 2 . 2 ) 


The result (2.2) has been obtained following a method by Brenner (1964) (cf. [Dorr fc 
Hardt (2014) concerning necessary corrections to the method), which is based on spherical 
harmonics expansions and yields the velocity and pressure fields around a slightly deformed 
sphere. To this end, the particle shape (given by the pair of fused spheres of radius a in 
our case) needs to be parameterised in spherical coordinates (r, 0, yi) according to 


r = rp{0,(f), 

and subsequently expanded in a power series 

rp{0, ip) = a 1 + (0, p) + {0, p) 


(2.3) 


(2.4) 


in terms of a small parameter e. Here, we choose e = cos 0, so that 2ea equals the distance 
between the centers of the fused spheres. Accordingly, if we assume the centers of the 
spheres to lie on the x-axis with the symmetry plane given by a; = 0, the particle shape is 
described by 


rp{0,p) = I 


1 + esin^ |cos(^| + e' 


iSm 


^ cos^ p — 1 


(2.5) 


from which the functions and in equation (2.4) can be read off. At the same 


time, the velocity and pressure fields, u and p, are written in the form 


u = 


eu 


( 1 ) 


■ £ 21 ^( 2 ) 


P — - 

The particle moves with the velocity U = Ue^ 
boundary conditions 

= u 

at the particle surface Sp, and 

m|v = 0 


+ 0{e^), and (2.6) 

-0{e^). (2.7) 

. Therefore, the flow field obeys the 

( 2 . 8 ) 


(2.9) 

on a spherical surface Eoo at r —>■ oo. While condition (2.9) is readily adapted to the 


perturbation expansion (2.6), condition (2.8) on the particle surface requires a Taylor 


series expansion for removal of the implicit dependence on the shape parameter e. To be 


precise, the boundary condition (2.8) in conjunction with the expanded particle shape (2.4) 
reads 


u{rp,0,p) = U, 


( 2 . 10 ) 
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which equals 


u''°'>{rp,0,(p) + eu^^\rp,0,ip) + e‘^u‘'^\rp,0,(p) +0{e^) = U. 


( 2 . 11 ) 


In equation (2.11), the argument depends on e, so that the velocity vectors {i G 
{0,1, 2}) are to be expanded in a Taylor series about r = a in powers of e, using 
equation (2.4). After performing the expansions, inserting the results into equation (2.11), 


and grouping of terms, we arrive at 


(o)=I7 


u 


= —ac, 


= —ac, 


,(i) 


,( 2 ) 


dr 

dr 


> at r = a. 


r 

2 L 


a(i) 


l2 


dp2 


— at 


a(i) 


dr 


( 2 . 12 ) 


Clearly, the zeroth order problem corresponds to a spherical particle moving with constant 
velocity U in an unbounded fluid, for which the velocity and pressure fields 


“-s'® ())'(( 




—|-3-)sin0ee, and (2.13) 
r a/ 


( 0 ) 3 COS0 

2 


(2.14) 


are well known ( Happel fc Brenner[|1983 ). The first-order problem has been solved by |Ddrr 
& Hardt (2014) in the particle’s rest frame. Because the frame of reference only affects the 
zeroth-order flow by a ddition or subtractio n of the velocity field I/e^, the first-order 
flow field considered by Dorr & Hardt (2014) may be directly used in the present study. 


The supplementary material contains the complete set of expressions requi red to calculate 
the velocity field u'A). Since is equal to the infinite series ( Brenner]|l964 ), 

the number of included terms needs to be limited in practical calculations. The values 
reported below as well as in the supplementary information correspond to /c ^ 20. With 
this choice, the (!l(e^)-contribution to the drag coefficient can be calculated to three 
significant digits, as will be shown in the following section. 


3. Applying the Lorentz Reciprocity Theorem 

According to the above discussion, we are in a position to compute the velocity 


fields (a, 0, (/?), zS {0,1,2}, on a spherical surface by means of equation ( |2.12[ ). In 
other words, we are faced with two flow problems involving a sphere of radius a in an 
unbounded fluid, see figure The first of these problems (figure [^a)), constructed by 
setting s to zero, consists of a sphere with a no-slip surface condition and translating 
with velocity I/e^. We shall denote the velocity and stress tensor fields belonging to this 
first problem by u and T, respectively. The second problem (figure [^b)), associated with 
the truncated perturbation expansion in e and denoted by u and T, is given by a sphere 
with a prescribed surface velocity field of 

OO 

u(o, 0, (f ) = ^ e®it^®^(a, 0,(p) = Ubz + eu^^\a, 0, ip) + e^it^^^(a, 0,(p) + 0{e^), (3.1) 

i=0 


according to equation (2.121. Since the solutions (u, T) and (u, T) correspond to the same 


flow geometry, they are related by the Lorentz reciprocal theorem for Stokes flow (Lorentz 
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Figure 1. Pair of flow problems connected by the Lorentz reciprocal theorem (3.2|: (a) a sphere 
moving along a fluid-fluid interface with vanishing viscosity ratio causes a flow field identical to 
the classical Stokes flow problem in a bulk fluid, denoted by u and i, if the contact angle 0 
equals 90°; (b) for 0 90°, the flow problem is equivalent to a pair of fused spheres moving 

through a bulk fluid; by means of a perturbation expansion, the boundary condition on the 
complex particle surface is projected onto a sphere; the flow field in this case is denoted by u 
and T. 


1896 Happel fc Brenne^|1983 ), 


(T-n)-udE= / (T-n)-udE, 


(3.2) 


where E comprises the particle surface Ep (outer normal vector n = —Er) and the 
surface Eco at infinity (outer normal vector n = e^.). The Lorentz theorem (3.2) has been 


exploited in a number of cases, for instance by Stone & Samuel (1996), Masoud & Stone 


(2014), and Schonecker & Hardt (2014). Its use in the present study is inspired by the 
argumentation of Stone & Samuel (1996). The contribution from the integral over Eq^ in 

(3.3) 


equation (3.2) vanishes because 

-1 


u 


-1 


r 


u 


T • Ej 


r and IIT • Er 


r ^ for r 


oo. 


The integral in the Lorentz theorem (3.2) thus reduces to an integral over the par ticle 
surface Ep. Recalling that T • n = —T • Er = 3fiiU/{2a)Ez ( Stone &: Samuel||l996 ) and 
using equation (3.1), the Lorentz theorem can be written as 

f , 


2a 


-Ue, 


UE,+eu<^^'>+0{e^) dE = / (T • n) • C/e^dE. 


(3.4) 


Since Ue^ is a constant vector and dE = 47ra^, equation (3.4) simplifies to 


c rr , 
GtthiU + 

2a 


eu 


(i)+eV2)+0(e3) 


dE = e. 


T • ndE. 


(3.5) 


= -Fr 
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Figure 2. Comparison of the models (2.2| and (3.7l with experime ntal and theoretical drag 
coefficient values taken from the literature. In the experiments by Petkov et al. (19951, an 
air-water interface with a very small viscosity ratio, ~ 0.02, was studied, 'i'he remaining 

curves are valid under the assumption = 0. 


The integral on the right-hand-side of equation (3.5) has been identified with the negative of 
the Stokes drag Fd on the particle because n = —Br- Using equati ons (2.12|), (| 2.13), (2.4), 
and (2.5) in conjunction with the velocity field according to Dorr & Hardt (2014), 


the integral occurring in equation (3.5) can be evaluated, yielding 


FD ■ = —GiTfiiUa 


16 


0{e^) 


(3.6) 


From the symmetry of the problem, it clearly follows that Fd = Fue^, so that the drag 


coefficient (1.2) for the original problem of a spherical particle diffusing along a fluid-fluid 


interface of zero viscosity ratio is given by 


/(0,o) = i 


1-1-cos 0 — 0.139 cos^ 0 -I- Ofcos^ 0) 

16 


Correspondingly, for the diffusion coefficient (|l.l|) it follows that 
D = 


16kT 


371/^10 [16 -I- 9 cos 0 -I- 2.224 cos^ 0 -I- O(cos^ 0)]' 


(3.7) 


(3.8) 


Note that the numerical coefficient 0.139 in equation (3.7) can be calculated to any desired 


number of significant digits, provided that a sufficient number of terms is considered in 


the spherical harmonics expansion by Brenner (1964) and Dorr & Hardt (2014). As stated 


above, the value 0.139... corresponds to a number of 20 terms. The underlying rational 
expression reads 765 368 413 099/5497 558138 880. 


4. Discussion 

In Hgure|^ we compare the result (3.7) and its first-order part (2.2) to experimental and 
theoretical drag coefficient values from the literature. The numerical solution to the full 
flow problem by Zabarankin (2007) may serve as a reference. Remarkably, equation (2.2) 
can be seen to be applicable at least down to a contact angle of 60°, where it still agrees 
well with the reference data. The inclusion of the cos^ 0-term, resulting in equation (3.7), 
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leads to a significantly better agreement between the asymptotic model and the numerical 
solution. Within a maximum error in the drag coefficient / of less than 0.02 (3%), the 
quadratic expression may even be applied to the full range of contact angles between 0 
and 90°. An approximate formula by Wang (2013), 


/w = 1 + 


sin(20) - 20 
^ ’ 


(4.1) 


based on the particle’s cross-sectional area seen by the respective fluid phases is also 
plotted in figure Apparently, this expression is at variance with the other four data 
sets shown in the figure. The region where 0 > 90° is excluded from figure owing to a 
lack of reference data. The drag coefficient is expected to approach zero when the contact 
angle approaches 180°, a behaviour which is not reproduced by the model equation (3.7). 


Therefore, the range of validity of the model for contact angles larger than 90° remains an 
open issue. In order to model the drag coefficient near 180°, either the above perturbation 
expansion would need to be extended to higher orders in e, or a reference body differing 
from a sphere (such as a disk or an ellipsoid) would need to be used. 


5. Conclusions 


With equations (3.7) and (3.8), we have developed explicit expressions for the drag 


and diffusion coefficients of a spherical particle attached to the interface between two 
immiscible fluids for the case of a small viscosity ratio between the two phases. The 
relations account for the dependence on the contact angle between the two fluids and the 
solid surface. Following from the assumption on the viscosity ratio, the drag and diffusion 
coefficients of a pair of fused spheres moving perpendicular to their line-of-centers has been 
found simultaneously. A comparison between the model and reference data has shown that 
the model can be applied to the entire range of contact angles below 90° with high accuracy. 
By applying the Lorentz reciprocal theorem to a geometric perturbation expansion, one 
order of approximation has been gained for the drag and diffusion coefficients as compared 


to the flow field. In other words, the second-order result (3.7) is based on the first-order 


velocity field, while calculating the non-trivial first-order result (2.2) only requires the 


well-known flow field around a sphere. The method can be applied to any particle shape 
resulting from a small geometric modification of another particle shape with a known 
flow field. 


Acknowledgments 

We gratefully acknowledge financial support by the German Research Foundation 
through Grant No. HA 2696/25-1. 


REFERENCES 

Ally, J. & Amirfazli, A. 2010 Magnetophoretic measurement of the drag force on partially 
immersed microparticles at air-liquid interfaces. Colloids and Surfaces A: Physiochem. Eng. 
Aspects 360, 120-128. 

Blawzdziewicz, j., Ekiel-Jezewska, M.L. & Wajnryb, E. 2010 Motion of a spherical particle 
near a planar fluid-fluid interface: The effect of surface incompressibility. J. Chem. Phys. 
133, 114702. 

Brenner, H. 1964 The Stokes resistance of a slightly deformed sphere. Chem. Eng. Sci. 19, 
519-539. 

Brenner, H. & Leal, L.G. 1978 A micromechanical derivation of Pick’s law for interfacial 
diffusion of surfactant molecules. J. Colloid Interface Sci. 65 (2), 191-209. 











Aaron Dorr and Steffen Hardt 


Chen, W. & Tong, P. 2008 Short-time self-diffusion of weakly charged silica spheres at aqueous 
interfaces. Eur. Phys. Lett. 84, 28003. 

CiCHOCKi, B., Ekiel-Jeewska, M.L., Nagele, G. & Wajnryb, E. 2004 Motion of spheres 
along a fluid-gas interface. J. Chem. Phys. 121, 2305. 

Danov, K., Aust, R., Durst, F. & Lange, U. 1995 Influence of the surface viscosity on the 
hydrodynamic resistance and surface diffusivity of a large Brownian particle. J. Colloid 
Interface Sci. 175, 36-45. 

Danov, K.D., Gurkov, T.D., Raszillier, H. & Durst, F. 1998 Stokes flow caused by the 
motion of a rigid sphere close to a viscous interface. Chem. Eng. Sci. 53 (19), 3413-3434. 

Dorr, A. & Hardt, S. 2014 Driven particles at fluid interfaces acting as capillary dipoles. 
arXiv: 1411.1183vl [cond-mat .soft]. 

Du, K., Liddle, J.A. & Berglund, A.J. 2012 Three-dimensional real-time tracking of nanopar¬ 
ticles at an oil-water interface. Langmuir 28, 9181-9188. 

Fischer, T.M., Dhar, P. & Heinig, P. 2006 The viscous drag of spheres and filaments moving 
in membranes or monolayers. J. Eluid Mech. 558, 451-475. 

Fuleord, G.R. & Blake, J.R. 1986 Force distribution along a slender body straddling an 
interface. J. Austral. Math. Soc. Ser. B 27, 295-315. 

Happel, j. & Brenner, H. 1983 Low Reynolds number hydrodynamics. The Hague: Martinus 
Nijhoff Publishers. 

Lorentz, H.A. 1896 Eene algemeene stelling omtrent de beweging eerier vloeistof met wrijving 
en eenige daamit afgeleide gevolgen. Zittingsverslag van de Koninklijke Akademie van 
Wetenschappen teAmsterdam 5, 168-175, translation into English: Journal of Engineering 
Mathematics 30: 19-24, 1996. 

Masoud, H. & Stone, H.A. 2014 A reciprocal theorem for marangoni propulsion. J. Eluid 
Mech. 741 R4, 1-7. 

O’Neill, M.E., Ranger, K.B. & Brenner, H. 1986 Slip at the surface of a translating-rotating 
sphere bisected by a free surface bounding a semiinfinite viscous fluid: Removal of the 
contactline singularity. Phys. Fluids 29, 913. 

Peng, Y., Chen, W., Fischer, Th.M., Weitz, D.A. & Tong, P. 2009 Short-time self-diffusion 
of nearly hard spheres at an oil-water interface. J. Fluid Mech. 618, 243-261. 

Petkov, J.T., Denkov, N.D., Danov, K.D., Velev, O.D., Aust, R. & Durst, F. 1995 
Measurement of the drag coefficient of spherical particles attached to fluid interfaces. J. 
Colloid Interface Sci. 172, 147-154. 

POZRIKIDIS, C. 2007 Particle motion near and inside an interface. J. Fluid Mech. 575, 333-357. 

Radoev, B., Nedjalkov, M. & D.iakovich, V. 1992 Brownian motion at liquid-gas interfaces. 
1. Diffusion coefficients of macroparticles at pure interfaces. Langmuir 8, 2962-2965. 

Schonecker, C. & Hardt, S. 2014 Electro-osmotic flow along superhydrophobic surfaces with 
embedded electrodes. Phys. Rev. E 89, 063005. 

Sriram, I., Walder, R. & Schwartz, D.K. 2012 Stokes-Einstein and desorption-mediated 
diffusion of protein molecules at the oil-water interface. Soft Matter 8, 6000. 

Stone, H.A. & Samuel, A.D.T. 1996 Propulsion of microorganisms by surface distortions. 
Phys. Rev. Lett. 77 (19), 4102-4104. 

Walder, R.B., Honciuc, A. & Schwartz, D.K. 2010 Phospholipid diffusion at the oil-water 
interface. J. Phys. Chem. B 114, 11484-11488. 

Wang, D. 2013 Tracer diffusion at water-oil interfaces studied by fluorescence correlation 
spectroscopy. PhD thesis, Universittsbibliothek Mainz. 

Wang, D., Yordanov, S., Paroor, H.M., Mukhopadhyay, A., Li, C.Y., Butt, H.-J. & 
Koynov, K. 2011 Probing diffusion of single nanoparticles at water-oil interfaces, small 
7 (24), 3502-3507. 

Zabarankin, M. 2007 Asymmetric three-dimensional Stokes flow about two fused equal spheres. 
Proc. R. Soc. A 463, 2329-2349. 



9 


Appendix A. First-order velocity field 

A.l. Series representation of the velocity field 
The first-order velocity field can be calculated according to the method by 
(1964). However, as pointed out by Dorr & Hardt (2014), corrections to Brenner’s method 
are necessary. Here, we provide all of th e expressions r equired to explicitly compute the 
velocity field which is of the form (jBrenner]|l964 ) 


Brenner 


u 


( 1 ) _ 


= E 


u 


( 1 ) 


(Al) 


k=0 


To facilitate the explicit evaluation of the velocity field and to simultaneously ensure a 
high degree of accuracy, we truncate the series after 20 summation terms and arrive at 


,(i) 


20 


,(i) 


The velocity fields can be expressed by 


(A 2) 


u 


where 


( 1 ) _ 


OO r 

E Vx 


n — 2 


2n{2n — l)/i 




I ^ Ji) 


( 1 ) _ {2n-l)fii 

kP—n—1 




k^P-n-l 


n+1 


(n + 2),AW+,yW 


2(n -I- 1) Vr 

( 1 ) _ ^ 
kX-n-l - 


(n + l)a 


(“) 


n+1 


^(1) 

k^n 


(A3) 


(A 4) 
(A 5) 
(A 6) 


(Brenner 19641. Through accommodating the velocity field (A3) to the boundary condition 


at the particle surface, the functions kXn \ kXn^'^ and kZn’ can be specified. One has 

(A 7) 


.TTT 


=0 for n < 1, 


y(i) = 

n 




2 ^ 2k + I 

m—Q 
k 

3 CkmNkm,cos{nnp) 


for n = fc — 1 


- S 


m—0 


2k+ 1 


lo 


(2 -I- fc)(fc — m -|- l)P^;^(cos0) for n = fc -I- 1 

for all other n, 
(AS) 


and 


^(1) _ J ^UJ2^^Qmckm,NkmPir'{cos 0)sm{m(p) for n = fc 
I 0 for n ^ k. 


k n 


(A 9) 


In equations (A7)-(A9), the P™ are associated Legendre polynomials. Expression (A8), 


developed by IDorr & Hardt (2014), differs from the corresponding result by Brenner 
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(1964), while equations (A7) and (A9| have been adopted from Brenner (1964) without 


modification. The expansion coefficients Ckm and are given by 


hk+l {k-m)\ -r 


m > 0 
if m = 0 


(A 10) 


and 


TT 277 


^km — 


J j sin0 |cos(/?| A^fcm-P™(cos0) cos(m(/?) sin0d0d(/3, (All) 


0 0 


respectively. Equation (All) contains information about the particle geometry via the 
term sind |cos (^|, which is equal to the first-order term in the expansion of the particle 
shape (equation (2.5) in the paper). The velocity field according to equation (A 1) is 
thus fully determined. For the sake of convenience, we explicitly display the coefficients Ckm 
for fc ^ 20 in appendix |A.2| Due to the symmetry of the particle shape with respect to 
the plane cc = 0, coefficients with m < 0 and/or m odd vanish. 


A.2. Coefficients c^r 


k=0: m=0 


k=2: m=0,2 

_ \/57r 

8 ’ 8 


k=4- m=0,2,4 

_ 3\/^ y/'oT^ VSStT 

64 ’ 32 ’ 64 


k=6: m=0,2,4,^ 

5vT^ 3V9l¥ 

1024 ’ 1024 ’ 1024 ’ 1024 


k=8: m=0,2,4,6,8 

SSCiTtv 3-y/S^ 3Vl3097r 3Vl21557r 

16384 ’ 4096 ’ 8192 ’ 4096 > 16384 


k=10: m=0,2,... ,10 

_147^ 49^/38^ 

131072 ’ 131072 ’ 65536 ’ 131072 ’ 131072 ’ 131072 


k=12: m=0,2,... ,12 

_3465^ 45V3603if 225^^1^^ 75^^^ 15Vl385677r 15y/ 15^^71^ 

1048576 ’ 524288 ’ 1048576 ’ 524288 ’ 1048576 ’ 524288 ’ 1048576 


k=14: m=0,2,... ,14 

14157V297r 1089V395857r 33^4688183571- 33^1278595571- 33^5881539377 

33554432 ’ 67108864 ’ 33554432 ’ 67108864 ’ 33554432 ’ 67108864 ’ 

_ 165-v/r^SmI^ 495-^/64 632 371- 

33554432 ’ 67108864 
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k=16: m=0,2,... ,16 

30673573311 2044979357r 3003746 1 897r 


1073741824 ’ 

268435456 > 

268435456 ’ 268435456 ’ 


1001^/"53im6l 
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